Withholding Potentials, Absence of Ghosts and Relationship 
between Minimal Dilatonic Gravity and f(R) Theories 
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We study the relation between Minimal Dilatonic Gravity (MDG) and f(R) theories of gravity and 
establish strict conditions for their global equivalence. Such equivalence takes place only for a certain 
class of cosmological potentials, dubbed here withholding potentials, since they prevent change of 
the sign of dilaton $. The withholding property ensures the attractive character of gravity, as 
well as absence of ghosts and a tachyon in the gravi-dilaton sector and yields certain asymptotic of 
the admissible functions f(R). Large classes of withholding cosmological potentials and functions 
f(R) are found and described in detail. It is shown that the popular choices of f(R) functions are 
not withholding ones. The particle content of the gravi-dilaton sector is found using perturbation 
theory around de Sitter vacuum of MDG. The graviton remains massless, since it obeys conformal 
invariant field equation in the de Sitter space-time. The R/Q term in the conformal invariant wave 
operator introduces a very small mass scale tur « 1.5 x 10~ 38 m e , m e being the mass of the electron. 
The mass of the dilaton is much larger: m 2; 10 29 m_R. Two new phenomena: scalaron waves and 
induction of gravitational waves by the scalaron field are discussed using the derived wave equations 
for scalaron and graviton. The MDG and f(R) theories are shown to predict physical deviations from 
GR. Seemingly, the MDG and f(R) theories, when globally equivalent, offer a unified description of 
dark energy and dark matter. 

PACS number(s): 04.50.Kd, 04.62. +v, 98.80.Jk 

I. INTRODUCTION 

Consider a scalar-tensor model of gravity with action 
for the gravi-dilaton sector [l|-[5[ [IH : 



Ag,* = Yk U 



j d A x^f\g~\(<S>R-2kU{<5>)). (LI) 



Without any relation with astrophysics and cosmology, 
this model was studied as early as in [l[ to justify the 
Fujii idea about the "fifth force" 0- There the term 
"dilaton" for the field <E> was introduced jilj]. Therefore, 
we call this model the minimal dilatonic gravity (MDG) 
0-0 • It corresponds to the Branse-Dicke theory with 
identically vanishing parameter uj. This property is also 
used as a name of the model in some publications. 

Here n — 8ttGn/c 2 is the Einstein constant, is 
the Newton gravitational constant, A is the cosmological 
constant, and <E> £ (0, oo) is the dilaton field. The values 
$ must be positive since a change of the sign of $ entails a 
change of the sign of the gravitational factor Gn/& and 
leads to antigravity instead of gravity. Such a change 
is physically unacceptable. Besides, the value $ = 
must be excluded since it leads to an infinite gravitational 
factor and makes the Cauchy problem in MDG not well 
posed [9j . The value $ = oo turns off the gravity and is 
also physically unacceptable. 
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The function <7(<£>) defines the cosmological potential 
which must be a positive single valued function of the 
dilaton field $ by astrophysical reasons[65j. Further on, 
we shall justify additional physical requirements for the 
cosmological potential [/($) as a necessary ingredient of 
a sound MDG model. 

The scalar field $ is introduced to consider a variable 
gravitational factor G($) = Gat/$ instead of the gravi- 
tational constant Gat. The cosmological potential f/($) 
is introduced to consider a variable cosmological factor 
instead of the cosmological constant A. In general rela- 
tivity (GR) with cosmological constant A we have $ = 1 
and (7(1) = 1 . Due to its specific physical meaning, 
the scalar field <J> must have quite unusual properties in 
comparison with the well-known classical spinless fields 
described in detail, for example, in [Tol [TT|. 

Some physical and astrophysical consequences of MDG 
are described in 0-H1, H|. In @ a theory of cosmological 
perturbations for general scalar-tensor theories, including 
MDG, was developed as a generalization of the approach 
of [8] . A recent attempt to develop a theory of cosmolog- 
ical perturbation using directly action (|I.1[) can be found 
in [l|. 

In particular, in [Hll it was shown that MDG can de- 
scribe simultaneously: 1) The inflation and the graceful 
exit to the present day accelerating de Sitter expansion 
of the Universe. 2) The reconstruction of the cosmologi- 
cal potential [/($), using the time evolution of the scalar 
factor a(t) in the Friedman-Robertson- Walker model. 3) 
The way to avoid any conflicts with the existing solar 
system and laboratory gravitational experiments using a 
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large enough mass of the dilaton field m$ > f 0~ 3 eV/c 2 . 
Masses above this one, even very large: m$ > Mpi anc k 
are not excluded by the existing experiments. 4) The 
time of inflation as a reciprocal quantity to the mass of 
dilaton m$. 

In the last decade the f(R) theories of gravity with the 
action 

A m) = d 4 xV\g~\f(R) (1.2) 

attracted much attention, since they seem to offer a pos- 
sible explanation of the observed accelerating expansion 
of the Universe without introducing "dark energy" and 
may be also an explanation of the observational problem 
of missing mass without introducing "dark matter", see 
fl3T - [45| and a large amount of references therein. For 
f(R) = R — 2 A we obtain once more GR with the cos- 
mological constant. 

Some of the popular choices of the function f(R) are: 

f . The Starobinsky 1980 model f(R) = R + 

^ 2 R 2 , n > 0. 

2. The Carroll et al. 2004 model [16]: f(R) = R - 
H 4 /R, n > 0. 

3. The Capozziello et al. 2006 model [20]: f(R) = 

ee (0,1). 

4. The Appleby et al. 2007 model [21]: f(R) = R + 
R c In (e-» + (1 - e~») e - R / R ") , fi, R c > 0. 

5. The Hu et al. 2007 model [13]: f(R) = R- 
fiR c (l + (R/R c r 2r ^ , n,R c >0. 

6. The Starobinsky 2007 model [H,[I1: f(R) = R- 

nR c (i - (i + {R/Rc) 2 ) n,R c >0. 

7. The Amendola-Li et al. 2007 model (Hi, [H]: 
f(R) = R- fiR c (R/R c y, »,R c >0,pe (0, 1). 

8. The Tsujikawa 2008 model [13 : f(R) = R - 
fxR c tanh(R/R c ), ^,i? c >0. 

9. The Appleby et al. 2009 model [29] f(R) = 
R + c(eln{^^)-R) + l(R/M) 2 , e = 

Rc 

b+ln(2ch(b)) ' 

10. The Gannouji et al. 2012 model (jO): f(R) =R- 
[iR c (l - (l + (R/Rcf^ H ^j+aR 2 , n,a,R c >0. 

In the above cases the function f(R) = R + Af(R) 
is constructed as a GR action plus corrections. Some of 
these functions f(R), being carefully chosen, are qualified 
as cosmologically viable. Their further modification and 
consistency with the available astrophysical data is still 
under investigation and may be considered as a promising 
approach to the modern astrophysical problems. 



Many authors think that using the Legendre transfor- 
mation [H, of the / (inaction JO) to the MDG- 
action (|I.1|) and vice versa they are able to prove the 
equivalence of these two theories of gravitation, both of 
which generalize GR. We show that a careful analysis 
leads to the opposite conclusion and refutes the claims 
that the f(R) theories are, in general, physically equiv- 
alent to MDG. From mathematical point of view these 
two models are, in general, equivalent only locally. The 
local equivalence is not enough to ensure identical physi- 
cal content and the same physical predictions of the two 
models. 

Therefore in the present paper we study a new prob- 
lem: The global equivalence between MDG and the f(R) 
theories. We derive additional requirements and define 
the class of the cosmological potentials E/($) which in- 
deed lead to a global equivalence and thus avoid some 
of the well-known problems, like physically unacceptable 
singularities, ghosts, etc., in the f(R) theories. Some of 
these requirements are known and dispersed in a large 
amount of the existing literature on the f(R) theories, 
where the problem of the global equivalence with MDG 
has been never discussed to the best of our knowledge. 
One of the goals of the present paper is to study this prob- 
lem, to collect all results at one place, to represent their 
strict derivation, and to describe all necessary and suffi- 
cient conditions for their usage, thus making transparent 
the mathematical structure and the physical content of 
the theory. 

An obvious inequality of the MDG and f(R) theories 
lies in the fact that we do not have a physical intuition 
on how to choose the function f(R). In contrast, our 
large experience in different kinds of physical theories, 
starting with classical mechanics, is very helpful in the 
construction of the cosmological potential C/($) 

In Section [TT] we present the derivation of different 
forms of the field equations of MDG and a discussion 
of their specific features and peculiarities. 

In Section IIIII we consider the specific Legendre trans- 
form between MDG and f(R) theories and derive the 
conditions for the global equivalence of these two mod- 
els. 

In Section llHBI we demonstrate the advantage of MDG 
introducing a new type of withholding dilatonic poten- 
tials V($) and cosmological potentials U(Q). A basic 
goal of this paper is to find large enough classes of such 
potentials that may be able to describe real problems 
(Section 1MB 2p . We also derive the corresponding form 
of a withholding function f(R). 

In Section ITVl we demonstrate the solution of the ghost 
problem in both the MDG and f(R) theories using the 
witholding potentials. 

The form of a withholding potential V(<p) in the Ein- 
stein frame is derived in Section HVBI 

Another result, presented here, is a correct formulation 
of the perturbation theory around de Sitter vacuum in 
MDG and clarification of the particle content of the gravi- 
dilaton sector and its basic physical properties (Section 
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IIV C[) . Here we use once more the withholding property 
of MDG. 

In the Concluding Remarks (Section[V| we outline the 
basic results of the present paper and some open physical 
problems. 

II. DERIVATION OF THE FIELD EQUATIONS 
OF MDG 

We derive the field equations of MDG using the stan- 
dard variational principle for action 1|) in some four- 
volume with a three dimensional boundary dV^ 1 . 

A. Variation of the action with respect to the 
dilaton field $ 

The variation of the action (|I.lj) with respect to the 
dilaton field $ is extremely simple and gives 

S^A g ^ = I d 4 x^\(R- 2AZ7,*($))<5$, (Ill) 
2k J 

[/$ being the derivative with respect to the field 
Hence, without any restrictions of the variation (5$ on 
the boundary dV 1 , using the variational principle one 
obtains the algebraic relation 

R = 2AU (II.2) 

A specific feature of MDG is the existence of such a 
relation, instead of a differential equation for the field <E>. 
According to (|II.2[) , the dilaton field $ has no space-time 



If the variations of the metric coefficients <5g Al/ and the 
variations of their derivatives S (<9„.g Aly ) on the surface 
dV 1 ^ obey the system of four restrictions: 

<5>5v x - V„$ 5g Xv + g Xv V„$ 6 (In \g\) = 
= 6 fag*) + d» (M^VWl)) + T aP x Sg^ + 
+g x "(d v (6]n\g\)-(d l ,$)6ln\g\) = 0, (II.4) 

then the standard variational principle yields the vacuum 
equations for gravy-dilaton sector in the form 

®G aP + AU($)g a0 + V a V^ - g a0 n$ = 0. (II.5) 



properties and evolution independent of the space-time 
scalar curvature R. In particular: 

i) In space-times with R = const due to relation (|II.2|) 
we have $ = const and gravitational factor G($) — 
const. This corresponds to our physical expectations for 
gravity, for example, in a flat space-time with R = 0. 

ii) In the case of validity of the cosmological princi- 
ple (CP) at very large scales (~ 10 2 ~ 10 3 Mpc) the 
3D space is homogeneous and isotropic, maybe even flat. 
Then all quantities must depend only on the cosmic time 
t. As a result of CP, we have R = R(t). Then (jlL2|) 
gives $ = $(£), and the gravitational factor automati- 
cally obeys the CP. 

In any other scalar-tensor model, or in any other field 
frame, instead of the algebraic equation (|II.2|) we will 
have a partial differential equation for the dilaton field 
<&(x), due to the existence of a kinetic term for $ in the 
corresponding action. As a result, the field $ would be 
a new scalar degree of freedom with an independent evo- 
lution. Thus, we would lose the above attractive prop- 
erties. These seem to correspond to the physical reality 
with good precision, when one interprets the dilaton field 
$(x) as a local strength scale of gravitational interaction. 
For example, it would be quite strange to have some dy- 
namics of the local intensity of gravitational interaction 
in an empty flat space-time. 

B. Variation of the action with respect to the 
metric field g a p 

The variation of the action 1[) with respect to the 
metric field g a p is much more complicated. Taking into 
account the identity (|A.5[) (see Appendix A) we obtain 



(II.3) 



The trace of the eqs. (|II.5[) (taking into account Eq. 
(III.2[0 . yields the following dynamical equation in vac- 
uum: 

□$ + AV>($) =0. (II.6) 

By construction it is written in terms of the dilaton field 
$. Here V>($) = § - 2U) = |$ 3 ^ (<S>- 2 U) is 

the derivative of the dilatonic potential. For convenience, 
we use the normalization V(l) = 0. Then 

2 /"* 

V($) = ■x / ($17,* - 2U) d$. (II.7) 
Ji 

Independently of pi.7[) . we suppose that the very point 
$ = 1 obeys the equation V|$($) = 0. As a result, we 



g A g ,t> = J cPx^f^G^ + AU($)g a0 + V Q V /3 $ - g af} U^j 6g a/3 
■ j> a?a x ^/\g~\ (^>5v x - V„$ Sg Xu + g Xv V <5 (In \g\) 



c 

2k 
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have U$(l) = 2 0]. 

The other independent field equations of the system 
(III.5P are defined by its traceless part [13] in the form 

<f>RP = -V^V^$. (II.8) 
If we add the standard action of the matter fields VP: 

A ma tt = ~f <Px^J\g\L matt (*, V*; Soj9 ), (II.9) 

which is based on the principle of the minimal interac- 
tion with gravity and contains only the metric, we obtain 
for the gravi-dilaton sector the dynamical system (in cos- 
mological units A = 1,k = 1,c=1 which will we use 
further on everywhere): 

□$ + \ ($£/,*($) -2f/($)) = \t, (11.10a) 
<S>R f l = -vlv^$ - T£, (II.10b) 

being the energy-momentum tensor of the matter. 

In this case, the dilaton $ does not interact directly 
with the matter and is a good candidate for the dark 
matter. Its interaction with the usual matter goes only 
trough the gravitational interaction. In the language of 
QFT this interaction is a process of second order with 
respect to the Newton constant, i.e. a very weak one. 

Now we are ready to define the de Sitter vacuum (dSV) 
as a state in which 

T£ = 0, * = $ = 1: V,*(1)=0, => 

Gt = Si (11.11) 

and the Einstein vacuum (EV) as a state in which 

Gi=0, $ = $ : tf,»(«o)=0, 

if + U(Q> )8i = 0. (11.12) 

A remarkable feature of the whole system (|II.2|) , 10|) 
is its invariance under the local transformations 

{U, T} h> {U + u{x),T -Au{x)} (11.13) 

with an arbitrary local field u(x) which is independent of 

There is one more circumstance that needs a comment. 
Actually, Eq. pi.lObp is the simplest form of the dynam- 
ical equation for the space-time scalar curvature i?, not 
for the dilaton <E>. Indeed, assuming 7^ one obtains 
from (1TTT2T) 

9 = f iR (R). (11.14) 

Here / R (R) is the first derivative of some new function 
f(R). Substituting this relation in Eq. pi.lObj) we obtain 
it in an explicit form 15a,|) . If in addition f RR 7^ 0, 
Eq. (|II.15a[) is equivalent to Eq. pi.l0a[) . Besides, under 



the same conditions we can exclude the dilaton $ from 
Eqs. (jll.lObj) . as well. Thus, we obtain the new system 

f >RR UR + f RRR (Vi?) 2 - i (Rf iH -V) = \ T, (IL15a) 

f, R Rl = -/. M Cv*fl- f, RR nV^RV?R-Tl (II.15b) 

The dynamical equations (|II.15j> of the problem look 
more complicated than the system (jll.lOp . The advan- 
tage of the form (|II.15|) is in the absence of the dilaton 
field $ in it. If one solves this system, the field <I> can be 
obtained using relation pi. 141) . The disadvantage of the 
simple form (|II.10|) is that it generates a wrong feeling 
about possible "independent dynamics" of the dilaton <E>. 
This is a little bit subtle circumstance. The character 
of Eq. pi,10a[) . which actually describes the dynamics 
of the scalar curvature R, can be explained in the con- 
text of the Legendre transform. A specific property of 
this transform is that the basic relations take the sim- 
plest possible form, if one uses a mixed representation in 
which half the variables are new and the other half is the 
old ones. 

The physical equivalence of the above two approaches 
to the field dynamics in the gravi-dilaton sector of MDG, 
based on Eqs. (|II.10|) . or Eqs. pi,15[) . needs a careful 
study. If such equivalence takes place, Eqs. pi.l0a|) and 
(I11.15a[) present two equivalent forms for description of 
the additional scalar degree of freedom which comes into 
being in both the MDG and f(R) models. It is frozen 
in GR where the scalar curvature does not own an in- 
dependent degree of freedom since R — T. Following 
paper [23], we call this new field degree and the related 
with it spinless particle the scalaron (For more details see 
Section ]IVC]). 



III. THE LEGENDRE TRANSFORM, THE 
RELATION BETWEEN MDG- ACTION AND 
/(-R)-ACTION, AND THE WITHHOLDING 
PROPERTY 

A. The Legendre transform between f(R) theories 
and MDG. 

Relations pL2]) and (|TlT4"f show that the MDG and 
f(R) theories are related via the Legendre transform, 
described, for example, in [4(| and [47} • Following the 
traditional notation in gravity, in our problems we ob- 
tain some particular form of the Legendre transform. We 
present here some of the specific details. 

For equivalence of the MDG and f(R) theories their 
actions p.lj) and (|I.2|) must give equivalent results under 
the corresponding variations. For this purpose, we need 
to satisfy the relation J(i?, $) = f(R) + 2C/($) - i?$ = 
d a v a {x) with some vector field v a (x) which does not 
dependent on R and $ [6q . It is enough to have 
^(i?,^) = f(R) + 2[/($) - i?$ = const. Then such a 
vector field v a (x) certainly exists. 
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In any case the conditions d<j,3 r (R, $) = and 
d R J(R, $) = produce relations pL2l and pLMl) , re- 
spectively, i.e. the function ^(R, $) generates the Legen- 
tre transform. Adding the condition|69| <J>) = 0, we 
obtain the transformation from t/($) to f(R) - pil.la[) . 
and the inverse transformation f(R) to - (Tffl.lbj) . 

in the following parametric form: 

/ = 2($[/,J$) -[/($)), 

i? = 2t/ s ($), $€ (0,00), (IILla) 

$ = / >H (ii), .Re (-00,00). (Ill.lb) 

The relation $,_r-R $ = 1 and Eqs. yield the 

new ones: 

I.rrU^ = ^, /,$ = 2$£/<j,$ = 
C^fi = \Rf : RR = \R$,R, 
UU R = f, R U* = \R®, (III.2) 

which will be used as a dictionary for translation of the 
results in MDG to f(R)-results and vice versa. 

All the above relations are mathematically correct and 
certainly physically acceptable if and only if we have con- 
vex functions: £/($) - in the interval $ G (0, 00), and 
f(R) - in the interval R € (—00, 00). Let us remind the 
definition: a twice differentiable function g(x) is convex 
in a given interval x € (x m i n ,x max ) if and only if its 
second derivative g"(x) ^ everywhere in this interval. 

Applying the Legendre transform to a function which 
is not convex, we will obtain a multi-valued new function. 
This means that the correspondence between MDG and 
f(R) theories will be only local. To what extend such a 
local equivalence may be acceptable from a physical point 
of view is a problem which we discuss in more details in 

Section niEll 

Hence, there is a one-to-one global correspondence be- 
tween the MDG and f(R) models if and only if the 
functions [/($) and f(R) are convex in their physical 
domains. Only under this condition the second of the 
equations Eqs. (|III.la|) and (jULlbj) can be solved unam- 
biguously, and the two models are physically equivalent. 



on the dilaton potential V(§). Then the infinite potential 
barriers at the ends of the physical domain (jIII.3l) will 
confine dynamically the dilaton ^>(x) inside this domain, 
if it is initially there [J, |5|. We call this new phenomenon 
the withholding property of the dilatonic potential V(&). 

For further use we need to derive the consequences of 
this condition for the cosmological potential t/($). For a 
given dilaton potential V(&) under normalization U (1) = 
1 the (IIII.5I) gives 

[/($) = -$ 2 / Q-Zv&d® + $ 2 . (III.5) 
2 Ji 

Let us choose a behavior of the dilatonic potential 
V(4>) ~ v $" (with some constant v > 0) to satisfy the 
condition pil.4[) at the ends of the interval $ g (0,oo). 
Then from Eq. (|III.5[) we obtain: 

a) For $ — > 0: n < 0, and 

[/($)~^rH-l (111.6) 

2 \n\ + 3 

b) For $ — > 00: n > 0, and 

{$ 2 for n e (0,3), 

|v$ 2 ln$ for n = 3, (III.7) 
l^a"*"- 1 for n > 3. 

As seen, in each case we have 

U(0) = U(oo) = +00, (III.8) 

but the increase of E/($) and V(§) at the ends of the 
physical domain (|III.3|) is, in general, not the same. 

An observational astrophysical fact: the cosmological 
term has a definite positive sign in the observable Uni- 
verse and leads to the additional requirement 

£7($) > for $ G (0,oo). (III.9) 

For potentials with properties (|III.8[) . (|III.9[) the convex 
condition reads 



B. The withholding property 

1. The withholding property in the interval $ € (0, 00) 

Looking at the dynamical Eq. (|II.6|I for dilaton $(2;) = 
$ (R(x))) we see that the only way to force the dilaton 
to stay all the time in its physical domain 

$€(0,oo), (III.3) 

preserving each value $ > attainable, is to impose the 
conditions 

V{0) = V{oo) = +00 (ILIA) 



> for $ e (0,oo). (IIL10) 

It ensures the uniqueness of the Einstein vacuum. We 
call such U ($) the withholding cosmological potentials. 

It is not hard to derive the qualitative behavior of ad- 
missible functions f(R) of general type created via the 
Legendre transform of the withholding cosmological po- 
tential [/($). Indeed, from Eqs. (hTI.6|) and (hTI.7|) we 
obtain the correspondence: $ —> R —> —00, / — > 
—00, and $ -> +00 <^> R -> +00, / -> +00. The 
translation of the properties (|III.3|) and (|III.9[) in the 
language of the f(R) models is f R > and f iRR > 
for all R € (—00,00). Hence, f(R) is for sure a strictly 
monotonically increasing and convex function. 



G 



For cosmological potentials U with asymptotic (|III.6[) 
and pil.7|) one obtains from pil.l[) : 



f(R) 



R^-oo 

M 



(IH +2) 



3-y In I 



l-i.- 



\R\ 



1 



n 



e (1/2,1), 



(in. li) 



f(R) 



R-^+oo 



r fi? 2 , for n G (0,3), 

9w$(i?) 2 ln$(i?), for n = 3 

( n _ 2 )(^3) 1/Al+ " 1 fA 

forn > 3, 



1/M+ (111.12) 



where in $(i?) = Lambert W (exp (yf^j)) we use the 
Lambert- W-function [Hj|, and ju+ = G (1/2, 1). 

The general form of the function f(R), which follows 
the asymptotic (|III.11[) and 12[) and is in addition a 
convex function, is shown in FigfTJ We dub such func- 
tions the withholding f(R) functions. 



Qualitative behaviour of a withholding function f(R) 



f 




/(«)-«! 



FIG. 1: Qualitative behavior of a typical withholding f(R) 
function 

The uniqueness of the physical de Sitter vacuum in 
such MDG is still not guaranteed. Indeed, from Eq. 
(|H.7|) we obtain 



= 3$[W(IIL13) 



These relations show that in the physical domain (|III.3|> 
the functions J7 1##s and V.^^^ have the same signs and 
zeros, but the functions {7 )4>4 , and do not own this 
property. Hence, V(&) may have several minima in the 
domain (|HL3)) [4], see Fig. |2f70j . 

Thus, we obtain a new kind of models with many 
locally stable dSV at the points $fe=o > 1, 
($o = 1) which are solutions of Eq. $fcf/ $ ($fe) - 



uv 




FIG. 2: Withholding potentials U (in blue) and V (in red) 
which have an unique EV and many dSV 



2U{$ k ) = and have U,„(* k ) > 2(7($ fe )/($ fc ) 2 > 



0. This means that the dilaton behaves around differ- 
ent minima like a field with different masses m\ k — 
(2/3) ($kU,**($k) ~ U,*($kj) > and excludes a scalar 
tachyon. In the language of the f(R) theories we 
have a series of dSV-curvature-values R k defined by Eq. 
Rkf, R(Rk)+f(Rk) = with /, R (R k )/R k > f, R R (R k ) > 
0. In each dSV state we have different values of the grav- 
itational factor Gfc = Gn/^u and of the cosmological 
constant A k = AU(®k) > A. 

These models do not need to introduce more new phys- 
ical structures and fields. An interesting questions is 
whether it is possible to construct a new mechanism 
based on dilatonic potentials with many dSV which is 
an alternative to the chameleon one and, in addition, de- 
scribes a drastic decreasing of cosmological constant we 
need for consistency with observations. 



UV 




FIG. 3: Withholding potentials U (in blue) and V (in red) 
which have unique EV and dSV 
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2. Examples of the f(R) models which are equivalent to 
MDG with unique dSV 

There is one more option: To postulate the unique- 
ness of the dSV [4|. In this case the function V r ($) 
with properties (|III.4[) will be convex: (<£>) > for 
each $ £ (0,oo). This ensures the stability of dSV and 
may exclude a scalar tachyon since m 2 = V ss ($) = 
(2/3) (E/ 44 (l) - 2) > 0, see Section [EVC] 

Besides, the function (2/3) ($[/^($) - = 
(/, R (R) - Rf, RR (R)) /3f, RR (R) > is strictly positive 
in the whole physical domain and this may avoid the 
presence of a scalar tachyon in the chameleon models, 
see reviews in (39l - l45| and the references therein. In ad- 
dition the condition /, rr(R) > entails the requirement 

f, R (R)/f, RR (R)>R, Re (-00,00). (111.14) 

For R = R = 4 this yields the condition for stability of 
dSV The stability of dSV in MDG is an obvious 

consequence of the withholding property. 

One can write down a simple example for a pair of 



By construction V(l) = 0, V R (1) = 0. The additional 
requirements 

w v+tV _ > 0, w »+>»- = 1 (HI- 18) 

ensure that f/(l) = 1 and V RR (1) = p~ 2 , as well as 
I 



/ = 



R = 



In all of these cases the necessary and sufficient con- 
dition for equivalence of the MDG and f(R) models are 
satisfied. The functions U(&) and f(R) are both con- 
vex and single-valued. The proper choice of the powers 



withholding potentials with unique dSV: 

Vm = t;P~ 2 + ^ - 2) , *e(0,oo), 

u@) = $2 + ^ 2 ( <I) -^) . ( IIL15 ) 

where p is a small parameter, related with the nonzero 
mass of the dilaton field[71] [4, 5], see FigJ3] 

The corresponding function f(R) can be written in the 
following parametric form[72|: 

/= ^(* 2 -| + 2 )+ 2 * 2 > 

R = 4^( $_ i) +4 *' ( IIL16 ) 

More general potentials with the same property ware 
described in [J] and can be used as a basis for a very 
general class of withholding potentials [73j [74j : 



(III.17a) 



(III. 17b) 

I 

that the convex condition U RR (&) > for 4> 6 (0, 00) is 
fulfilled. 



The corresponding withholding function f(R) is de- 
fined in the parametric form: 



(III.19a) 



(III.19b) 

I 

v± which enter into the sums can make the withholding 
barriers at the ends of the physical domain (|III.3[) impen- 
etrable not only in the classical dynamics of fields under 
consideration, but also at quantum level. 



_ 2 ^ W V+ , V _ $- v --l\ 



2<i> 2 + l V Wv ^- f(y+-^W + - 1 -^ K + 2)$— 1 + $ 2 

p 2 ^ U + + V. 



4$ 



p 2 ^ u+ + v. 



v+ — 3 V- + 3 

SI/+-2 off. I-,, 1 1\tf>— 2 



[v+ - l)^^- 2 - 2$ (u- + 1)$" 



2$ 



v + — 3 



V- 



v± > 0, 

$ e (0,00). 



s 



Much more strong withholding barriers at the ends of 
the physical domain (|III.3|I are produced, for example, 
by the potential 

a > 0, $ G (0,oo). (111.20) 

It yields a withholding function f(R) defined in a para- 
metric form 



3. Direct Construction of withholding functions f(R). 

The second relation in Eqs. pil.laP can be solved ex- 
plicitly in rare cases. Even if this is possible, as a rule, 
the result is complicated and not very useful for fur- 
ther usage. Taking into account the asymptotics (|III.11I) , 
12|) . and the convex condition we can construct the 
withholding functions f(R) in a more direct way. 

Consider the functions 



/ = 2 $ 2 + JL (( 2a ($ 2 - l/$ 2 ) - 1) e^*- 1 ^ 2 + l) , 

R = 4$ + -4 $4 ~ 1 e Q (^ 1 / $ ) 2 , $G(0,oo). (111.21) 
4p 2 a $ 3 

Its qualitative behavior is the same as the one shown in 

Figru 

From an analytical point of view, when in the formu- 
lae piI.17[) . piI.19[) only a finite number of terms en- 
ters, at the ends of the physical domain (|III.3|I we will 
have poles[75] of the corresponding functions, or at most 
log singularities - in cosmological potential t/(<E>) , and a 
more singular term related with the Lambert-W function 
- in f(R). Otherwise, the ends of the physical domain 
are essential singular points, as in the example pil.201) . 

mrm . 



F(x; fi, a) = (x 2 + 1) m/2 / (1 + cxp(-2aa;)) 

with the parameters /i, a and construct a five parametric 
family of f(R) functions 

f(R;Hi,H2;a,Ri,Ra) = 
= F(R/R V ,l/fi u a) - F(-R/R 2 ;fi 2 ,a) (111.22) 

with the parameters fix.2 G (1/2,1), and a,Rx,R 2 G 
(0, +00). By construction these functions have the 
asymptotic behavior 1 1|) . (|III.12|) and satisfy the 
Starobinsky condition /(0; fii,fi2j 01, Rx, R2) = 0. The 
next step is to construct the linear combinations, analo- 
gous to Eq. piI.17|) : 



}f(R;fi 1 ,^ 2 ;a,R 1 ,R 2 ), (III. 23a) 

E W{^i,M2,a,iii,ii2} = !j W{^ ljfl2>ajRuR2 } > 0. (III. 23b) 

{fi 1 ,li 2 ,a,R 1 ,R 2 }eCxD 

I 



Unfortunately, the problem of the convex property of 
these functions is not trivial. 

Let us call the convex domain (CxD) the domain of 
the parameters fj,i } ua, a, Ri, R 2 for which the functions 
pll.221) are convex |76l|. The CxD turns to be not a sim- 
ple set of points in the five dimensional parameter-space. 
If we perform the summation in Eq. pil.23a|) only over 
the values {fix, fi 2 , a, Rx, R 2 } G CxD, then the condi- 
tions pil.23b[) ensure the convex property of the func- 
tions f(R) and these functions are withholding. 

Thus, we see one more advantage of the MDG in com- 
parison with the f(R) theories: In MDG the convex do- 
main CxD has a simple form described in the relations 
piI.17p by the inequalities v± > 0. In contrast, in the 
f(R) theories CxD is a quite complicated set of points in 
the space of the parameters. 



4- The withholding property in a finite interval 
$ G ($1, $2) C (0, 00) 

Another possibility is to withhold the dilaton field in a 
finite interval $ G [$1, $2] where < $1 < 1 < <J> 2 < 00. 
We shall not discuss in detail the corresponding class 
of withholding potentials U and V. The withholding 
potentials with unique dsV of this class can be writ- 
ten in the form V := v w Vl<l , 2 V Vl< u 2 i and U := 
2i/i,i/a w »l,^ Uuiw ^1 ,2 6 (0 , 1), with coefficients w VuVa 
which obey relations pil,18[) . The dilatonic potentials 
V Vl>V2 have the form 

V Vl , V2 ($) = ^x (111.24) 
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0.5 1 1.5 2 



FIG. 4: Withholding potentials U (in blue) and V (in red) in 
the finite domain ($i, $2) = (0.5, 2) for 1/1,2 = 1/2 

where 6^,^ = ^i/(l — $1) — ^2/(^2 — 1) and c VuV2 = 

(i^i/(l — *&i) 2 + ^2/(^2 — I) 2 ) \ see an example in 
FigJU All these potentials satisfy the requirements 
V{1) = 0, V .(1) = 0, and V^(l) = 1/p 2 > and 
contain enough free parameters for fitting real problems. 
The potential E/($) can be found using Eq. (|III.5|I . 

It can be shown that the qualitative behavior of the 
corresponding function f(R) is the same as shown in 
Figffl but now the asymptotics are 

f(R) ~ ^R, f(R) ~ $ 2 i2, (111.25) 

R— > — OO ii— >+oo 
< $1 < $2 < OO, 

and corresponds to the cases with n± — > 1 in QUI. 111) . 
(IIlLT2|l . 



5. /(S,) models which are not equivalent to the MDG. 

1. For the Starobinsky 1980 model one has t/($) = 

^(<f-i) 2 , $ e (-00,00). 

2. For the Carroll et al. 2004 model the potential 
[/(<&) = ±/i 2 \/$ - 1, $ € (1, 00) is a double- valued 
function. 

3. For the Capozziello et al. 2006 model [/($) = 
§ ($/(l + e)) 1+1/e , $6(0,00). 

4. For the Appleby et al. 2007 model £/ = ^ ($ In $+ 
(1 - $)ln(l - $) - $ln(exp(/i) - 1) + ln(l - 
exp(-/*))), $e [0,1]. 

5. The Hu et al. 2007 model HU - the func- 
tion f(R) is never convex and the potentials 
[/(<&) qre multi-valued. For example, £/($) = 

2^ (V7*± a/1 = $) 2 : * e (-oo,l) - for n = 1. 
The case n = 2 is shown in the Figj5l 

6. For the Starobinsky 2007 model the parametric rep- 
resentation of the multi- valued potential V($>) can 
be found in [50j. The multi- valued potential ?/($) 
is shown in FigJS] 




0.5 1 1.5 2 



FIG. 5: The multi- valued potential E/($) of the Hu et al. 2007 
model (n=2) 




0.5 1 1.5 2 



FIG. 6: The multi-valued potential [/($) of the Starobinsky 
2007 model 

7. In the Amendola-Li et al. 2007 model we ob- 
tain [/($) = -p)( w ,/(l-$))( 1 -P)/P,$ g 
(-oo,l). 

8. For the Tsujikawa 2008 model we obtain the two- 
valued potential, defined in a finite interval: 

[/($) = -(R c /2) ( x /T _ $arccosh(^/i/(l - $))± 

±vV($-i + m)) , *e(i-M,i). 

9. For Appleby et al. 2009 model the potential 
£/($),$ e (-00, 00) is shown in FigfT] 

10. The Gannouji et al. 2012 model represents some 
improvement of the Starobinsky 2007 model. It has 
a CxD for large values of the parameter a. The two 
possible cases are shown in Figs J5] and [9j 

The models 1, 5, 6, 7, 8 (for /i > 1), 9 and 10 admit neg- 
ative values of dilaton $ and, in general, contain ghosts. 
The relations (|IV.2bj) have not meaning for all values of 
$, admissible in these models. The requirements (|III.4I) . 
(IIII.8I) are also not satisfied. 

In the models 6 and 9 one can choose a special initial 
condition $0 f° r the dilaton in a tiny domain, which (at 
classical level) ensures that the solution $(i) > is pos- 
itive for all time instants. Depending on the choice of 
the initial conditions the existing singularities in model 
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FIG. 7: The potential [/($) of the Appleby et al. 2009 model 
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FIG. 8: The single-valued potential [/($) of the Gannouji et 
al. 2012 model for a = 0.6 



ture (in cosmological scales) when the model maybe must 
be changed [78]. 

The functions f(R) of the models 2, 5, 6, 8, 10 (for 
small a) are not convex and yield multi-valued cosmo- 
logical potentials {/($), thus being not acceptable. 

The models 3 and 4 are not withholding because do 
not prevent the inadmissible value $ = 0. 

As a result, all the above models 1-10 are not compat- 
ible with MDG. 

In the existing literature we can find also other choices 
of function f(R). For example, terms ~ R~ n (hiR/fi ) , 
~ exp(-W?)-l, or~ (exp(bR)-l)Hex£(bR)+exp(bRo)) 
can be included, see the reviews [35|, |42j and the refer- 
ences therein. In the known to us variety of such func- 
tions we were not able to find a f(R) model which is 
consistent with the withholding property of MDG. 

IV. PRESENCE OR ABSENCE OF GHOSTS 

In the present paper we call "a ghost" any physical 
field with a wrong sign of the kinetic term in the corre- 
sponding lagrangian[79j. The problem of ghosts in the 
so called theories of gravity with higher derivatives has 
been studied in many articles, see for example, the arti- 
cles dHH, EH El, the review articles and the 
references therein. For exclusion of ghosts in the Branse- 
Dicke theory with nonzero parameter ui and massless, or 
very light dilaton, see [^|[80j. 



6 can be not reachable for certain class of classical so- 
lutions. This shows that the proposed in the present 
paper requirement for a global equivalence of MDG and 
f(R) theories is acceptable but not obligate at classical 
level [73|. However, such a fine tuning can not prevent 
evolution to negative values of $, or to a negative cos- 
mological factor, if one takes into account the quantum 
tunneling. 

Another possible scenario to avoid the negative val- 
ues of $ only at the present epoch is to choose a proper 
positive value of $0 and to tune the parameters of the 
model in a way which ensures positive $(t) > for long 
enough time, i.e. to postpone the disaster for distant fu- 




0.5 1 1.5 2 





FIG. 9: The multi- valued potential U(&) of the Gannouji et 
al. 2012 model for a = 0.08 



A. The conditions for the ghost elimination in 
MDG 

It is not hard to derive the necessary and sufficient 
conditions for the absence of ghosts in the MDG. 

1 . There is no ghost-problem for dilaton field since in 
the action (|I.1[) a kinetic term for $ does not exists 
at all. 

2. The problem with the presence of ghosts in the 
metric field g a p is nontrivial. According to Sec- 
tion 94 in [6|, in the gauge go a = 0, a = 1,2,3; 
^g = const we can choose the coordinates in the 
Riemannian space-time of the theory to be Carte- 
sian at an arbitrary fixed point. Then, at this point 
9a/3 = <? Q/3 = —Sap for a, (3 = 1, 2, 3, and in the ac- 
tion (|I.1[) we will have under the invariant integral 
only terms of the form 

i$.g 00 (9 .9^) 2 . (IV.l) 

The kinetic term for metric field in p.l[) consists 
of terms pV.l[) multiplied by the positive constant 
c/2k. Since g 00 > 0, the necessary and suffi- 
cient condition for the absence of ghosts in MDG 

is mra . 
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The above results show that the conditions (jIII.3|) and 
pil.4l) , which define the withholding properties, ensure 
in a dynamical way also the absence of ghosts in MDG. 
Hence, the withholding property is a critical condition 
for the physical consistency of MDG. 



B. Presence of ghosts in f(R) theories of gravity 
and their relation with MDG 

The widespread opinion is that despite the presence 
of the second derivatives of the metric g a p in the f{R)- 
action (|I.2[) (see the last equation in (|A.2|) and the def- 
initions (|A.1[) '), there is no problem with ghosts in this 
modification of GR (see, for example, the recent paper 
[3lj). As described in some detail in [32], the standard 
reason for this belief is based on the following two steps: 

1) Legendre transformation of the /(i?)-action (11.21) to 
the MDG-action (JUTJ) . 

2) Conformal transformation of the MDG-action 1[) 
to the Einstein-frame-action: 



Aef = 



-ld*xV\~9\R- 



(IV.2a) 



2ld A x^/\f\(^d al pd^-V(4,)), 



g a i3 = ^g a p, <p = 

V(<t>) = e - 2 VV~^U 




(IV.2b) 
(IV.2c) 



which resembles the Hilbert-Einstein-action of GR, plus 
an additional scalar field ip with a normal kinetic term 
and self-interaction V(ip). Obviously, ghosts are not 
present in the Einstein- frame-action pV.2a|) . However, 
the action pV.2a[) is produced by the two relations 
(|IV.2b[) which both lose meaning outside the physical 
domain pil.3[) . Hence, the accurate analysis shows that 
the f{R) theories of gravity are not automatically free 
of ghosts. For this to be true certain strong additional 
requirements equivalent to the corresponding constraints 
on the MDG are needed. 

Some authors mention that for the absence of ghosts 
in the f(R) theories one needs the conditions 



f, R (R)>0 and f, RR (R)>0. 



(IV.3) 



However, the natural question: What kind of physical 
requirements can ensure the fulfilment of these conditions 
and thus indeed can exclude ghosts, remains without a 
definite answer. We saw that these conditions do not 
guarantee the withholding property of the function f{R), 
as well as that the withholding property is critical for 
exclusion of ghosts. For this purpose, relations pV.3[) 
are not enough, and the function f{R) must have proper 
asymptotic behavior, say, described by Eqs. and 
(HTTT2]) . 

As a result, in Einstein-frame-action (|IV.2I) the poten- 
tial V(ip) must have asymptotic 



and for ip — > +00: 

V(<f)r 

A simple example which corresponds to the cosmolog- 

(IV.6) 



const for n G (0, 3), 

ip for n = 3, (IV.5) 

e (n-3)y/2/3 V for n > 3 _ 



ical potential (|III.15|) is 
3 



V(ip) = l + ^p- 2 [l-e-^V^ 
For cosmological potentials we obtain 




„(|n|+3)V 2 / 3 M 



for 



ip ■ 



(IV.4) 



Unfortunately, the requirement for withholding prop- 
erty is largely ignored in the current literature on the 
f(R) theories, see the Section fill B 51 

Another cosmologically valuable example is V(ip) = 
Vo exp(— V3ip) (in our normalization). It was considered 
in [19| . There it was shown that this very simple potential 
offers a good fitting of the observational date in a large 
domain, showing some small deviations from standard 
ACDM model in the CMBR angular power spectrum for 
small I. Actually, in this domain both models do not fit 
very well the existing data. The situation is not clear due 
to the large observational errors, but for the exponential 
potential under consideration the deviations seem to be 
in a wrong direction. We would like to draw a special 
attention to this potential, because it is very close to the 
withholding property: In this case the dilatonic potential 
= const x $™, n = 3(1 - l/y/2) € (0, 1) does not 
exclude the point $ = with infinite gravitational factor, 
but excludes $ = 00. The cosmological potential U (<&) = 
const x <E>™ _1 leads to an infinite A-factor for $ = and 
to a zero limit for A-factor when $ — > 00. One may hope 
that a small change of this simple potential - just to obey 
the withholding property, may improve the fitting of the 
known observational data. 



C. Particle content of gravi-dilaton sector in MDG 
and its physics 

The particle content and the ghost problem of the f(R) 
theories was discussed as early as in [14|, see also the 

papered). 

In the present Section we consider the perturbations 
around dSV in MDG. Here the propagators have a more 
elegant form. The physical content of the model is more 
transparent and simpler than in the f(R) theories which 
are only locally equivalent to MDG, since a lot of possible 
cases, considered in [14(, are automatically excluded in 
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MDG. In the case of a global equivalence the physical 
content of the two models is the same. 

Consider small perturbation of the metric g^ v = + 
/i M „ around the solution g^ v for dSV ([Il.lljl . The linear 
approximation for perturbations of Eqs. (|II.10[) above 
dSV yields the particle content of the gravi-dilaton sector 
of MDG. The needed basic geometrical relations are given 
in the Appendix B. Indeed: 

1 . Let us introduce the scalaron £ as excitation of the 
dilaton above dSV in the following way: 



> „ j XD $ Dh-2h 



6 to 2 + 2C//3$ 



(IV.8) 



Then it obeys the scalaron wave equation with a 
source: 



□C + \R + m 2 C = 8T/3, & 



(ST/3), 



□ + ±i? + m 2 



(IV.9) 



where the mass of the scalaron, defined as 



iR for to 2 > iR 

6 — 6 



(IV.10) 



is real, and ST is the perturbation of the trace of 
energy-momentum tensor of matter. 

2. For the traceless part of the metric perturbation 
h a p (which carries the graviton[82j) we obtain the 
wave equation: 



(□ + iR) h aP 

h a 8 = = 



i(- 



ST 



a/3 



□ ■ 



-R $ 



vT^c) ■ (iv.ii) 



ST a p being the perturbation of the traceless part of 
energy-momentum tensor of matter. We see that: 

• The propagators of the scalaron (spin s = 0) 
and graviton (spin s = 2)[83j: 



n. 



□ ■ 



\R + m 2 



n 



s=2 



□ ■ 
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(IV. 12a) 
(IV. 12b) 



look like propagators of a standard classical 
massive field (scalaron, s = 0) and massless 
field (graviton, s = 2) in de Sitter space-time, 
correspondingly j84|- In the graviton case the 
additional term ^R ensures the conformal in- 
variance of the Eqs. pV.lip . and the zero 
mass of the graviton(85| . 

If m 2 > ^R, the scalaron never becomes 
a tachyon. This inequality can be fulfilled 



as a result of the withholding property since 
to 2 = V, $$(!>) > for $ being a point of dila- 
tonic potential minimum. We shall see that 
this condition is met with a large stock, since 
from another reason physically admissible are 
the values m 2 ^> hR. 

• MDG predicts the existence of free spreading 
scalaron waves. Since the mass of the dilaton 
is not zero, their velocity is smaller than the 
velocity of light. 

• The electromagnetic field and other sources 
with T = are not able to create scalaron. 

• As a result of the withholding property the 
graviton never becomes ghosts since $ > 0. 

• The scalaron is a source of gravitational waves. 
This phenomenon is natural for a theory (like 
GR) in which every physical field curves the 
space-time and thus becomes a source of grav- 
itational waves. 

• If one tries to express the scalaron contribu- 
tion in the source of gravitational waves di- 
rectly via ST/3, one will be forced to use the 
operator 



1 



□ 



1 



:V Q V £ 



1 



□ + 



iR- 



(IV. 13) 



i 



iR- 



After some manipulations one obtains for the 
second term in the brackets, enclosed between 
the two propagators, a representation which 
resembles a presence of a ghost. Such inter- 
pretation is certainly not correct since this op- 
erator is not a propagator but just a source- 
operator. Its physical meaning is clear: the 
matter creates scalaron waves via ST/3 ^ 
0. These scalaron waves induce gravitational 
waves. 

Thus, the particle content of the gravi-dilaton sector 
of MDG and the corresponding basic physics in it is clar- 
ified. 



V. CONCLUDING REMARKS 

In the present paper, we have studied carefully the re- 
lation between Minimal Dilatonic Gravity (MDG) and 
the f(R) theories of gravity. These models seem to of- 
fer a unified description of dark energy and dark matter. 
Our final result is that the two generalizations of GR 
are globally equivalent under certain additional assump- 
tions. Their physical equivalence takes place only for 
certain class of cosmological potentials introduced here 
and dubbed withholding potentials since they prevent a 
change of sign of dilaton $, thus preserving the attractive 
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character of gravity. It is shown that the same assump- 
tions ensure the absence of ghosts and tachyon in the 
gravi-dilaton sector of MDG. 

We studied in detail large classes of withholding cos- 
mological potentials and functions f(R) exploring their 
asymptotic properties and using the convex condition. It 
is shown that the popular choices of the f(R) functions do 
not own the withholding property and, as a rule, do not 
exclude ghosts, being only locally equivalent to MDG. 
Some of them suffer also from other shortcomings, dis- 
cussed in the text. 

We develop proper perturbation theory around de Sit- 
ter vacuum of MDG. The background space-time is de 
Sitter one. In this framework we defined the scalaron 
and graviton as a field excitation above the de Sitter vac- 
uum of MDG and justified their basic physical properties 
deriving the corresponding wave equations with sources. 

According to Eq. (|IV. 1 1[) . in MDG graviton is a mass- 
less particle in the de Sitter space-time. In the cor- 
responding conformal-invariant wave operator □ + ^R 
presents a constant term which resembles a mass term in 
a flat space-time. It corresponds to an extremely small 
mass, even in comparison with the electron ones m e : 

m R = (h/c)^jR/6= (h/c) a/2A/3 « 1.5 x lCT 38 m e .(V.l) 

This value is consistent with the known astrophysical 
data about the mass of the graviton (62|. Such term 
is needed to have a right quasiclassical limit of the the- 
ory, as well as to have a unique vacuum (up to a unitary 
equivalence] in the Fock space in QFT in the de Sitter 
space-time [54M30] . 

The mass of the scalaron to in Eqs. pV.9l) , (|IV.10jl may 
be arbitrary large. The only known restriction, needed 
to avoid conflicts with the solar system and laboratory 
experiments was pointed out in the Introduction. It gives 

to > 2 x l(T 9 TO e « 10 29 m R . (V.2) 

Thus, the second condition in Eq. (|IV.10[) is certainly 
fulfilled in realistic MDG models. As seen, in their gravi- 
dilaton sector we have two very different mass scales. 

We recover two new phenomena in MDG: 

1 . Scalaron waves freely spreading in vacuum with ve- 
locity much smaller than the velocity of light be- 
cause of relation (IV. 21) . This is consistent with the 
expected low velocity of dark matter particles and 
supports the hypothesis that the scalaron may be 
a good candidate for dark matter. An important 
problem is to find experimental methods and/or ob- 
servational tools for a direct registration of scalaron 
waves predicted by MDG. 

Our consideration shows that from a geometrical 
point of view the scalaron waves are freely spread- 
ing in vacuum perturbations of space-time scalar 
curvature R (See Eq. (|IV.8[) .). Thus, their exis- 
tence is an essential deviation from GR. Here the 
scalar curvature is rigidly related with the trace of 



the energy-momentum tensor of matter and van- 
ishes identically in vacuum. 

2. Induction of gravitational waves by the scalaron, 
according to Eq. (IIV. 1 1|) . Estimates of the mag- 
nitude of this effect related with astrophysical ob- 
jects and with the early Universe are needed to see 
whether this prediction of MDG may be of interest 
for the future more sensitive detectors of gravita- 
tional waves. 

According to the withholding condition, the dilatonic 
potential has at least one minimum. Hence, in MDG 
there exists at least one dSV. If there are several of them, 
as shown in Fig. [5J then the withholding condition en- 
sures the ordering < $o < $i < ■ ■ • < oo. 

Our physical normalization requires for the lowest min- 
imum $o = l and C/($o) = 1 — a natural choice for the 
case of unique dSV in MDG [3, Hj]. If there exist many 
dSV, this normalization corresponds to the hypothesis 
that the present state of the Universe is close to the low- 
est one. 

For different dSV, numbered by k = 0, 1, ... , the grav- 
itational factor has different values Gk — Gm/^u- The 
different values of the cosmological factor are Afe = AUk- 
Hence, as a result of the withholding property, the lowest 
dSV is a state with minimal value of the cosmological fac- 
tor and with maximal value of the gravitational factor. 
The scalaron has also different masses to^ in the vicinity 
of different dSV. 

This interesting new physical situation needs a more 
careful analysis. The model does not need introduction of 
additional physical structures and fields. An interesting 
open question is whether it is possible to construct a 
new mechanism based on dilatonic potentials with many 
dSV, which is an alternative to the chameleon one, and, 
in addition, describes a drastic decreasing of the initial 
cosmological constant which we need for consistency with 
observations. 

A basic open problems remains a careful choice of the 
withholding MDG model which is able to fit known ob- 
servational data in cosmology and astrophysics. 
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Appendix A: Derivation of the MDG-field-equations 

In this article, we use the following notation written in 
the sign conventions of the article [23| : 

r Q 0, 7 = 2 (dagpy + dpg ai — d 1 g a p) , 
V = 3 7 T a/3 , M , T Q = = 8 a (in v1s|) • (A.l) 

R a p,i 5 = —2 (9[ a r /3 ] 7 ' 5 + r[ Q ,| a | i5 r^] 7 ") , R a p = R^a,^ 1 , 
R = g^R^u = g*" (r M1 /r 5 ~ r M /r„ a ') + 

+^=d x (3"% - <nv, A )) . (A.2) 

The variation of the Hilbert-Einstein action uses the 
well-known relations 




G a /3 = Ra/3 — T<japR- (A.3) 

The explicit form of the vector 5v x is not needed in GR. 
The direct calculations give 

5v x = g Xv 8 g T^ - g aP 5 g Y af} x = 
= d v g Xv + T a0 x 8g^ + TJg x » + 2g x »5 g r„. (A.4) 

Then, after some more algebra we obtain 

V\g\Vx<S>Sv x = yf\g~\(v a V $ - 9aP U^j5g ap - 
-du (v^V A $ (6g x » - 2g Xv 8 (in (A.5) 

Appendix B: Geometrical perturbations of the de 
Sitter space-time 

We use the general formalism for perturbations of the 
metric g^ v — g^ v + h ul/ on a curved background space- 
time with metric g^ [6(. We have to make some changes 



in the basic formulas because of our conventions (See 
Appendix A). The perturbation = g^ v — g^ = 5g^ v 
is a small quantity of the first order. The sign s» denotes 
equalities valid in the linear approximation. The bar sign 
denotes the phone-space-time-quantities. 

The tensor of the affine deformation is 

r Q/ 3 7 « H ap i = | (W a h} + Vphl - VX/j) . (B.l) 
Then 

R a p-y S ~ -2 V[ ct i/ /3 ] 7 ' 5 , (B.2) 

and imposing the transversal gauge condition 

V>£ - \V a h = 0, h = h% (B.3) 

from Eqs. (jB.2[) and (|B.3[) one obtains the general form 
of the variations of the Ricci tensor and scalar 

5R a0 « iD^ + ^ij^-F^, (B.4a) 
SR ks lUh-h^R^. (B.4b) 

Taking into account the definition of the dSV in MDG 
it is easy to obtain the following dSV-values: 

$: V*($) = 0, ^($) = m 2 >0, (B.5a) 
U = [/($), U„ = 2(7/$, 

= 2U/$ 2 + §m 2 /$, (B.5b) 
R = 4, R a/ 3 — g a /3, 

Ra/3jS = — 3 (3a 7 5/3<5 — gaSijfij) ■ (B.5c) 

The corresponding values for /, f tR , and f tRR can be 
obtained using Eqs. 1 1) . 

As a result of Eqs. (|B.5I) we obtain the following final 
form of the perturbations above dSV-space-time: 

6 Rap ~ \ Oha/3 + f h a p + | Uhg a p, (B.6a) 
SR a p | □ h a/ 3 + | h a p — S^Rap^j + h a p, (B.6b) 

S(Rap) « |D/l^ + i^ = i (n+^)Aa/S,(B.6c) 

&R w | Dft — /i. (B.6d) 

Note that the used barred quantities represent the cor- 
responding values for any possible dSV. To simplify the 
notation, here we do not number explicitly the different 
dSV. 
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models of gravity (see, for example [52T | and the refer- 
ences therein) in the f(R) theories we are able to exclude 
ghosts under conditions (|IV.3f) . despite the presence of 
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